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It is shown, in D = 2 + 1 dimensions, that by merely imposing non-abelian gauge invariance on the temporal 
gauge ground state wavefunctional of an abelian gauge theory, a confining state is obtained. 



Confinement is a property of the vacuum of quantized non- 
abelian gauge theories. In the Hamiltonian formulation of 
quantized Yang-Mills theory, in D = d + 1 dimensions and 
temporal gauge, the vacuum is the ground state solution [A] 
of the time-independent Schrodinger equation 
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Physical states in temporal gauge, in SU(2) gauge theory, are 
required to satisfy the Gauss Law constraint 



(5^+^*4)^ = 



(3) 



which is the condition for invariance under infinitesimal gauge 
transformations. In strong-coupling lattice gauge theory, the 
ground state can be solved for systematically, order-by-order 
in powers of 1/g 4 1 1] 
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where the R n [U] are gauge-invariant expressions formed from 
Wilson loops on the lattice. If the sum in the exponent is 
truncated at order n = M, the resulting state solves the time- 
independent Schrodinger equation to M-th order in 1 /g 4 , but 
satisfies the physical state condition exactly. 

Similarly, we may guess that at weak couplings there is also 
a weak-coupling expansion for the ground state. This would 
have the form (in continuum notation, with lattice regulariza- 
tion implicit) 



V [A] -exp 
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where the Q n [A] are gauge invariant, and the approximate so- 
lution, with the sum truncated to some order in g 2 , satisfies 
the Schrodinger equation to that order. It is awkward to build 
the Q„[A] out of Wilson loops at weak couplings; the strong- 
coupling series suggests that such a construction would in- 
volve summation over infinite sets of loops in the continuum 
limit. Instead, we conjecture that that the building blocks of 
the Q„[A] are the covariant derivatives, and their associated 



covariant Green's functions. If that is so, then the choice of 
<2o is essentially unique. We know that in the g = limit, the 
ground state wavefunctional is just 
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Then a state which solves the Schrodinger equation to zeroth 
order in the coupling g, but satisfies the Gauss Law constraint 
exactly, is obtained by simply replacing ordinary derivatives 
by covariant derivatives; i.e. 
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where is the covariant derivative in the adjoint represen- 
tation of the gauge group. Field strengths are of course com- 
mutators of covariant derivatives, so this expression can be 
regarded as a functional of covariant derivatives only. 

The non-confining ground-state solution of the free field 
theory, eq. l|6}, is well-known. Many years ago, it was sug- 
gested by one of us |2|] that, at large distance scales, the effec- 
tive Yang-Mills vacuum state is simply 1 
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This vacuum state has the property of dimensional reduction; 
i.e. the computation of a space-like loop in d + 1 dimensions 
looks just like the computation of a loop in a li-dimensional 



of the 



Euclidean Yang-Mills theory. If both the vacua x 

(2) 

3 + 1 dimensional theory, and "PX of the 2+1 dimensional 
theory have this form, then the computation of a large planar 
Wilson loop in 3+1 dimensions reduces, in two steps, to a 
computation in two-dimensional Yang-Mills theory, i.e. 
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See also Halpern 01 and Mansfield Q. 
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In D = 2 Euclidean dimensions the Wilson loop can be cal- 
culated analytically, and we know that there is an area-law 
falloff. 

In Hamiltonian lattice gauge theory it is possible to calcu- 
late the vacuum state systematically in powers of 1/g 4 Q|, 
and the result agrees with the above conjecture for the contin- 
uum theory. Further support for this form of the ground state 
wavefunctional was obtained via lattice Monte Carlo simula- 
tions ||5|]. In an interesting paper from 1998, Karabali, Kim, 
and Nair 1 6] came up with a strong-coupling expansion for the 
vacuum state of the continuum Yang-Mills theory in D = 2 + 1 
dimensions (see also Leigh et al. |7|]). This state, to leading 
order in 1 /g 4 , also agrees with eq. (|8}. 

If we take eq. @ as the zeroth-order approxi mation to the 
weak-coupling ground state, and if the kernel 1 /y— 2? 2 would 
have a finite range, then the effective ground state at scales 
much larger than this range would have the dimensional re- 
duction form l|8}. In order to make this statement a little 
more precise in D = 2 + 1 dimensions, consider the decom- 
position of the field strength at fixed time, in the x — y plane, 
into gauge-covariant and gauge-invariant factors: 

Ff 2 ( X ) = mm 



♦"(*) = -Mr . /« 



12 



F? 2 (x)F{ 2 (pc) (10) 



In SU(2) gauge theory § a (x) is a unit 3-vector in color space, 
which can be made to point in any direction, at any posi- 
tion x, by a suitable choice of gauge. In the x — y plane, a 
gauge-invariant distinction between high and low-frequency 
field-strength components can be made in terms of the Fourier 
transform f(k) of f(x). For some reference momentum scale 
O, define 
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We will also write ^0 = e with 
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Now, in the spirit of self-consistent field approaches, let us 
approximate the gauge-invariant expression (j)(l/ \J — © 2 )(|) by 
its expectation value 



S(x-y) = U a (x) 
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Then f 8 (x) decouples from f (x), and the dependence of 
the wavefunctional on the long range, "low-momentum" com- 



ponent F^ 2 w ' a (x) of the field strength is given by the dimen- 
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sional reduction form d8}. The fluctuations of F, 
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then essentially local, since in d = 2 dimensions there is no 
Bianchi identity among field strengths which would induce a 
correlation. On the other hand, with a momentum cutoff at 
\k\ = O there is an intrinsic short-distance cutoff, on the order 
/ = tz/o, on the spatial variation of f low {x). Therefore, the 
dimensional reduction form for F™' only implies that the 
correlation length, extracted, e.g., from the correlator 



(Tr[(F<n*))lTr[(F/7 (y)f\) (Tt[(F{^ (x))^ (15) 

is less than /. This upper bound on the correlation length is re- 
duced as o is increased, until, at o w k/Ik, the Fourier modes 
S(k) at \k\ < a cease to be constant, and f low {x) does not fluc- 
tuate independently in regions separated by distances I < Ik- 
Thus the range Ik of the vacuum kernel is also an estimate of 
the correlation length, which is the inverse of the mass gap. 
Denoting ttik = 1% , we would then have niK ~ niQ+, where 
wo+ is the mass of the + glueball. 

At first sight, however, there would seem to be no reason 
that 1 /V— D 2 should have a finite range. 2 But in fact, for 
the covariant Laplacian in two Euclidean dimensions, there is 
good reason to believe that the corresponding Green's func- 
tion is finite range for almost any stochastic background A M 
field, no matter how weak. The relevant phenomenon is 
known as weak localization lIlOll . Consider a non-relativistic 
particle moving in a stochastic potential in one or two space 
dimensions. It is known that the eigenstates of the correspond- 
ing Hamiltonian H are all localized; i.e. they fall off expo- 
nentially outside some finite region, no matter how weak the 
stochastic potential may be. In that case, the inverse opera- 
tor has a finite range. Weak localization is expected to 
occur also in cases where the stochasticity lies in the kinetic, 
rather than the potential term, and this seems to be relevant 
to the case at hand, i.e. the two dimensional covariant Lapla- 
cian. If all of the ei genm odes of —2) 2 are localized, then both 
1 / (— © 2 ) and 1 /V— T> 2 are necessarily finite range. In fact, 
the range is finite even if there is only an interval of localized 
eigenmodes at the low end of the spectrum, with the bulk of 
the states non-localized lHllll2ll . 

The range of the vacuum kernel 1 j\J — D 2 is fixed by self- 
consistency. In two dimensions there is no Bianchi identity, 



If the covariant vacuum kernel 1 / V— ® 2 has a finite range 
Ik, then S(x — y) has the same range, and its low-momentum 
components are approximately constant, i.e. S(k) w y, for 
|£| <C k/Ik- If we also choose a <§; k/Ik, we find 



2 For this reason it was suggested by Samuel, in ref. 1 8] that one could in- 
troduce a variational parameter m 2 into the kernel, i.e. 1 / \f— D 2 +m 2 , as a 
way of introducing a finite range. The value of m 2 would be determined, in 
principle, by minimizing the energy expectation value of the trial vacuum 
state. 
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so a correlation between field strengths at different points in 
space, at equal times, can only arise from the non-local char- 
acter of the vacuum wavefunctional. As argued above, field- 
strength correlations are negligible on a scale at which the vac- 
uum wavefunctional has the form (|8j, and the wavefunctional 
(0 takes this form on a scale which is roughly equal to the 
range of the vacuum kernel. This means that the correlation 
length of field strengths in a typical vacuum fluctuation (i.e. 
the inverse of the + gl ueball mass), should approximate the 
range of the kernel 1 / \/—D 2 , but the range of the kernel, in 
turn, depends on the vacuum fluctuation. 

The first question is whether, in a vacuum fluctuation typi- 
cal of massless free field theory, having infinite range correla- 
tions among field strengths, the kernel can also have an infinite 
range. If the kernel has a finite range in such a background, 
then massless free field fluctuations are inconsistent; they do 
not arise from the behavior of the kernel that they imply. We 
check this behavior numerically in the following way: It is 
required to generate stochastically a gauge field background 
that would arise in a free theory, from a probability distribu- 
tion |1q e [A\ | 2 , where, in momentum space, 



L=32 



o 




FIG. 1 : The modulus of the vacuum kernel 1 / \J — © 2 in a "free field" 
background, but with various non-zero values of the Yang-Mills cou- 
pling g in the covariant Laplacian. The upper solid line is the ker- 
nel G(R) = v3/(2ti/?) corresponding to A"(x) = (or, equivalently, 
g = 0). Data at high g tends to fall near the lower dashed line. 
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The superscript T means "transverse". In this initial inves- 
tigation the transversality restriction is ignored, and the con- 
tinuum is replaced by a lattice. Let — V 2 be the (ordinary) lat- 
tice Laplacian, with eigenstates {(p„(x)} and eigenvalues {/„}, 
and let {.?„,,} be a set of random numbers taken from a nor- 
mal distribution with unit variance. Then a typical free-field 
vacuum fluctuation, deriving from the probability distribution 
I^Q^fA] | 2 (with the transversality condition dropped) is 



Combining color and space indices, the covariant lattice 
Laplacian can be regarded as a sparse 3L 2 x 3L 2 matrix. Stan- 
dard numerical routines are used to take the inverse square 
root of this matrix, to arrive at 
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To calculate the range, we form the gauge-invariant combina- 
tion 



n=2 I, 



(17) 



This construction is repeated three times, for each of the three 
color indices, to form the field A" (x). 

Although no coupling constant is used to generate the free 
A-fields, there is still a Yang-Mills coupling constant present 
in the definition of the covariant derivative. In lattice regular- 
ization, this coupling constant relates the A-fields to the SU(2) 
link variables 



Ui(x) =exp 



(18) 



From the link variables in the fundamental representation, we 
then obtain link variables A U in the adjoint representation 



A Uf{x) = \Tr[<y a Ui(x)o b u}(x)} 



(19) 



The covariant lattice Laplacian, in the adjoint representation, 
is then 



k=l 

(20) 



G(x,y) = JY,G ah (x,y)G^(y,x) (22) 

i.b 

and average G(x,y) over all points x,y on the lattice with fixed 
\x — y\ = R. Denote the result G(R). A linear fit to log[G(R)] 
vs. R determines the inverse m# of the range Ik of the kernel. 

The results, for couplings from g = 0.1 to g = 2.0 on a 
32 2 lattice are shown in Fig. [2 which plots G(R) vs. Sona 
semilog scale. The solid line is the result for g = in the con- 
tinuum, which is simply G(R) = ^/3/(2%R). Even at g = 0.1, 
where the average link variable (jTrfC/]) = 0.997 is very close 
to unity, there is a noticeable deviation from 1/R behavior, and 
at all the larger g values the 32 2 lattice is big enough to clearly 
see an exponential falloff in G(R). 3 

This data indicates that 1 / \/—D 2 is finite range, even in a 
massless free-field background. This implies that free fields, 
and probably any sort of non-confining vacuum fluctuations in 
2+1 dimensions, are inconsistent. The kernel is finite range 



3 We note that the data at g > 1.0 also approaches a limiting line, approxi- 
mately at the position of the lower line in Fig.Q This simply means that, in 
the way we have generated the link variables, taking g — > °° does not result 
in a random link distribution equal to the Haar measure. 
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FIG. 2: A comparison, at various lattice couplings (3 in Z) = 2 + 1 
dimensions, of the inverse range m% of the vacuum kernel, and the 
mass of + glueball reported in ref. 1 9]. Errorbars are comparable to 
symbol size. 

in general, dimensional reduction is obtained, and the simple 
vacuum wavefunctional is confining. 

Given that the approximate ground state wavefunctional 
•Pn of eq. is confining, the next question is how closely 
this state approximates the true vacuum state of SU(2) Yang- 
Mills theory. We have already argued that in the state 0, 
rriK ~ ttiQ+ . If is a good approximation to the true Yang- 
Mills vacuum state x V , Q Ue , then in the true vacuum state niK and 
Wq+ should still be approximately equal. To test this approxi- 
mation, we calculate the range of the kernel 1 / V— ® 2 in a 2D 
time-slice of a thermalized lattice, generated by Monte Carlo 
simulation of lattice Yang-Mills theory in D=3 Euclidean di- 
mensions. This procedure makes use of the fact that if 0[A] is 
an observable defined at fixed time t , then 

(O) = ^jDAO[A]e- s 

= (% rue \0\% rue ) (23) 

Comparison of mjc with the known lattice results for the mass 
of the + glueball is one measure of how well the zeroth- 
order state 1*0 [A] approximates the true ground state, since 
if *J>'™ e = ^Po, then niK ~ mn+. The procedure for calculating 
nig in the true vacuum is as follows: Link variables are gen- 
erated by standard lattice Monte Carlo of D = 3 dimensional 
lattice gauge theory with a Wilson action, and the covariant 
Laplacian is evaluated in a time-slice. From this one calcu- 
lates 1 / V— £> 2 and G(R), as described above. Then \ogG(R) 
is fit to a straight line in the range R G \k, k]; the kernel mass 
m.K is just — 1 times the slope. The procedure is repeated for a 
number (> 10) of independent thermalized lattice configura- 
tions, and the results are averaged. 

The values of mk obtained in this way at a number of dif- 
ferent lattice couplings p are shown in Fig.|2] where we have 
also displayed the + glueball masses that have been com- 
puted via lattice Monte Carlo by Meyer and Teper in ref. |9]. 
For the evaluations of niK we have used lattice extensions 



L = 16,24,32,40 at p = 6,9, 12, 18 respectively. With the ex- 
ception of p = 18, where Meyer and Teper used a slightly 
larger lattice (L = 50), these are the same lattice volumes used 
in ref. |9]. In the range of p studied, trig and m + appear to 
converge as p increases. This is the behavior one would ex- 
pect, if the true vacuum and the zeroth-order approximation 

to the ground state approach one another in the weak cou- 
pling, P — > °° limit. 

Despite this numerical success, it cannot be true that the 
zeroth-order state and the true ground state agree in all essen- 
tial ways. Since the zeroth-order state has the dimensional 
reduction form at large scales, the dependence of the string 
tension on the group representation of the quark-antiquark 
sources is the same as in D = 2 Euclidean dimensions, i.e. 
Casimir scaling. While this dependence is correct at interme- 
diate distance scales, and also asymptotically in the N — * °° 
limit, it just cannot be true asymptotically at finite N, where 
the asymptotic string tension depends only on the A^-ality of 
the quark sources. The simple zeroth-order state appears to be 
missing the center vortex (or domain) structure, which has to 
dominate the vacuum state at sufficiently large scales 1 13, 
So dimensional reduction is not the whole story at finite N, 
and higher-order, 1/N 2 suppressed terms in the vacuum wave- 
functional must come into play at sufficiently large distance 
scales. At present we have no systematic method for deter- 
mining these higher-order terms. 

It should be noted that there have been a number of ef- 
forts to determine the Yang-Mills vacuum wavefunctional be- 
yond weak-coupling perturbation theory. These are mainly 
variational in character; see in particular fl5l [Tr3[. A non- 
variational approach is followed by Leigh et al. |7] in then- 
recent work on the Yang-Mills ground state in 2 + 1 dimen- 
sions. There are some strong similarities between that work 
and the present paper. Leigh et al. arrive, by various manipu- 
lations, at a ground state which is the exponential of a bilinear 
expression, involving a finite range kernel connecting gauge- 
invariant variables related to the field strength. Due to the 
finite range of the kernel, the vacuum state derived in ref. |7] 
should have the dimensional reduction form when viewed 
at large scales. This, of course, is very reminiscent of in 
eq. 0. However, it is hard to compare these vacuum states 
directly. For one thing, the two states are expressed in terms 
of different variables; ref. |7] makes use of the Karabali-Nair 
variables in 2+1 dimensions. For another, the kernel of ref. 

01 is gauge-invariant, rather than gauge-covariant, and has a 
finite range even for a vanishing background field. That is not 
at all like the behavior of the kernel 1 /%/— £> 2 . Still, there 
are enough superficial similarities between these two states to 
make one suspect that they are related. 

We would like to conclude with some remarks about 3+1 
dimensions. The evidence presented above suggests that con- 
finement in 2+1 dimensions is simply a consequence of the 
physical state constraint (i.e. gauge invariance), plus the free- 
field limit of the vacuum wavefunctional. The key point 
is th at gau ge invariance introduces a gaug e covariant kernel 
1 / \J — D 1 , and this kernel, unlike 1 / \J — V 2 , has a finite range 
in almost any stochastic background, possibly due to the weak 
localization property in two Euclidean dimensions. The inter- 
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esting question is whether the gauge-covariant kernel still has 
a finite range, in a free-field background at fixed time, in the 
3+1 dimensional theory. It is known that the spectrum of — D 1 
has an interval of localized states in D = 4 Euclidean dimen- 
sions, in a confining background 11 llll7ll . But it is not known 
whether this operator has localized states in a free-field back- 
ground in D = 3 dimensions, or else whether the covariant 
kernel at D = 3 is finite range for some other reason. If the 
kernel is, in fact, finite range, then on large scales I^Pol 2 ef- 
fectively reduces to the probability distribution of a three Eu- 
clidean dimensional Yang-Mills theory, and the conjectured 
two-step dimensional reduction of eq. (|9jl, resulting in an area 
law falloff for large planar Wilson loops, would hold. The 
matrix operations would have to be carried out on 3L 3 x 3L 3 
matrices, rather than the 3L 2 x 3L 2 size that appears in 2 + 1 



dimensions. This additional computation cost, while restric- 
tive, is probably not prohibitive. 



Acknowledgments 

We thank Dmitri Diakonov, Robert Leigh, Adam Szczepa- 
niak, and Valentin Zakharov for helpful discussions. This re- 
search is supported in part by the U.S. Department of Energy 
under Grant No. DE-FG03-92ER4071 1 (J.G.), the Slovak Sci- 
ence and Technology Assistance Agency under Contract No. 
APVT-5 1-005704 (S.O.), and the Grant Agency for Science, 
Project VEGA No. 2/6068/2006 (S.O.). 



[1] J. Greensite, Nucl. Phys. B 166, 113 (1980). 

[2] J. Greensite, Nucl. Phys. B 158, 469 (1979). 

[3] M. B. Halpern, Phys. Rev. D19, 517 (1979). 

[4] P. Mansfield, Nucl. Phys. B418, 113 (1994). 

[5] J. Greensite and J. Iwasaki, Phys. Lett. B 223, 207 (1989). 

[6] D. Karabali, C. Kim and V. P. Nair, Phys. Lett. B 434, 103 

(1998) |arXiv:hep-th/9804132|. 

[7] R. G. Leigh, D. Minic and A. Yelnikov, arXiv:hep-th/0604060 
[8] S. Samuel, Phys. Rev. D 55, 4189 (1997) 

|arXiv:hep-ph/9604405|. 
[9] H. Meyer and M. Tepcr. Nucl. Phys. B 668, 111 (2003) 

|arXiv:hep-lat/0306019|. 
[10] B. Kramer and A. MacKinnon, Rep. Prog. Phys. 56, 1469 

(1993); 

P. Lee and T. Ramakrishnan, Rev. Mod. Phys. 57, 287 (1985). 



[11] J. Greensite, S. Olejmk, M. Polikarpov, S. Syritsyn 

and V. Zakharov, P hys. Rev. D 71, 114507 (2005) 

|arXiv:hep-lat/0504008|. 
[12] A. McKane and M. Stone, Ann. Phys. (NY.) 131, 36 (1981). 
[13] J. Greensite, Prog. Part. Nucl. Phys. 51, 1 (2003) [arXiv: 

hep-lat/0301023 ]. 
[14] J. Greensite, K. Langfeld, H. Reinhardt, S. Olejmk, and T. Tok, 

arXiv:hep-lat/0609050 
[15] I. I. Kogan and A. Ko vner, Phys. Rev. D 52, 3719 (1995) 

|arXiv:hep-th/9408081 |. 
[16] A. P. Szczepaniak and E. S. Swanson, Phys. Rev. D 65, 025012 

(2002) |arXiv:hep-ph/0107078|. ^ 
[17] J. Greensite, A. V. Kovalenko, S. Olejmk, M. I. Polikarpov, 

S. N. Syritsyn and V. I. Zakharov, arXiv:hep-lat/0606008 



